BEC in a star-comb graph 
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We investigate the properties of free bosons hopping on a star-comb network, discussing the 
single-particle spectrum and the main thermodynamic equilibrium properties: Bose-Einstein critical 
temperature, fraction of condensate, and spatial boson distribution. We find an enhancement of the 
critical temperature with respect to other inhomogeneous networks. 



Although it is well known that free bosons hopping on 
translationally invariant networks cannot undergo Bose- 
Einstein condensation at finite temperature if the space 
dimension d is less or equal to two, very recent studies 
QI S B 3 hint to the exciting possibility that the net- 
work topology may act as a catalyst for inducing a finite 
temperature spatial Bose-Einstein condensation even if 
d < 2. This is indeed possible if one resorts to a suitable 
discrete inhomogeneous ambient space on which bosons 
are defined. As an example of this situation we shall 
investigate the properties of non-interacting bosons hop- 
ping on a star-comb shaped network, evidencing that also 
in this case one may have a macroscopic occupation of 
the ground-state at low temperatures. Ultracold bosons 
in optical networks represent a physical system where 
this setup could be experimentally implemented. 

The star-comb graph is a bundled graph which can be 
obtained by grafting a star graph to each site of a lin- 
ear chain, called backbone (see Fig. [1]). The sites of the 
graph can be naturally labelled introducing three integer 
indices (x, y, z) with x,y,z € Z, where x = 1, ... ,p la- 
bels the different arms on each star (excluding the links 
connecting the different stars), y — 0, . . . , L represents 
the distance from the backbone and z = 1, . . . , N s labels 
the different stars. In the following we shall assume pe- 
riodic boundary conditions (x, y, 1) = (x, y, N s + 1). The 
Hamiltonian describing non-interacting bosons hopping 
on a star-comb graph can be written as: 
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In Eq. |T|). t is the hopping parameter while 
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is the creation (annihilation) operator for bosons; 

is the number operator at site 
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(x,y,z). The filling, i.e., the average number of parti- 
cles per site, is defined as / = Nt/N s (pL + 1), where 
Nt is the total number of bosons and N s (pL + 1) is the 
number of sites. The adjacency matrix of a star-comb 
graph is given by: 
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The eigenvalue equation —t^2-Aijip(j) = Eijj(i), with 






FIG. 1: The star-comb graph. 



the previous labelling of sites, reads: 

^ [(8y>,y-l + S y ', y +l) (1 - Sy,o)8 x>x >5 ZjZ > 
x',y',z' 

+ 3y,0 <V,l<5z,z' + Sy,0 Sy',oSx,x' {&z',z-l +5z',z + l)] 

■^{x',y',z') = Eip(x,y,z) 
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By exploiting the translation invariance in the direction 
of the backbone, a Fourier transform in the variable z 
reduces Eq. ([3]) to a 1-dimensional eigenvalue problem. 
If one defines: -^(x^y^) = '}2 z e lkz il)(x,y 1 z) 1 with k = 
2im/N s and n — 1,...,N S , the eigenvalues equation ^ 
becomes: 

— t [($y',y-l + $y',y+l) (1 ~ ^y,o)^x : x' + <5y,0 5y',l 
x' , y' 

+ 2 cos(k )S y<0 5y' t0 S x ,x']'>P(x',y') = Eip(x,y) 

(4) 

where ^(x^^) = S(k — ko)ip(x,y), with fco = 2im/N s , 
n = . . . N s . Equation ([4]) can be regarded as the equiv- 
alent problem of a quantum particle hopping on a sin- 
gle star and interacting with a potential in the center, 
V(ko) = — 2£cos(fco). The eigenvalues and eigenvectors 
of Eq. (|4]) are determined by requiring that the free par- 
ticle solutions on the arms are solutions also in y = 0. 
Indeed, for y ^ 0, the wave function satisfying Eq. (g]) 
is: 
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Requiring that - on each arm - tp(x,y) is a solution of 
the eigenvalue equation at y = L introduces restrictions 
on the parameters h, A x , B x , and on the eigenvalues E. 
In other words, it amounts to require that A x , B x and h 
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should satisfy the p equations 

= -2t cos(h) (A x e ihL + B x e" lhL ) ( 6 ) 

x = 1,- • • ,p. 

The eigenstates defined in Eq. ([5]) should also satisfy 
(p — 1) matching conditions in the center of each star. 
Thus, one has 

A x + B x = A x+1 + B x+1 x = l,---,p-l. (7) 

Using Eqs. ©, the condition in one of the centers gives 
one more equation 



5^ (A x e ih + B x e- ih ) + 2 cos(fc ) (A x + B x ) 
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-2icos(ft) (A x > + B x >) . 



with cc = 1, • • • , p. Equations ((6]) and ([Sj> may be grouped 
in a homogeneous linear system of 2p equations which al- 
lows to fix the 2p parameters A x and B x . Upon denoting 
with Q the (2p x 2p) matrix whose elements are the co- 
efficients of the linear system given by Eqs. © and (jHJ), 
and requiring that 

detg = 9(/i,L)(l-e 2l/l ( L+1 )) p - 1 ■ {(p-2)cos(ft) 

, (9) 
-pcot [h(L + 1)] sin (ft) + 2cos(fc )j 

is equal to zero, the solution is unique. In Eq. © 
|0(ft, L)\ = 1 for any value of ft. 

One immediately sees that the values of ft for which 
ft = nn/{L + 1) (with n — 1,2,...,L) provide a set of 
L(p — l)-fold degenerate eigenstates of Eq. © for each 
value of fco (in total N s sets) . In addition, the solutions 
of the transcendental equation 

{p - 2)cos(ft) -pcot [h(L + 1)] sin (ft) + 2cos(fc ) = 

(10) 

provide the values of h associated to non-degenerate 
eigenstates. Equation (flU)) can be solved numerically and 
yields a set of N S (L — 1) non-degenerate eigenvalues cor- 
responding to different values of fco and to values of h 
which - in the thermodynamic limit - are equally spaced 
and separated by a distance tt/(L+ 1). As a result, each 
set of delocalized states is formed by (pL — 1) states cor- 
responding to energies ranging between — 2t and +2t. If 
the potential in the center U(fco) is equal to zero, one 
recovers the same condition of the single star [4|. 

Since the total number of states should equal N s (pL + 
1), there are still 2N S localized states in the spectrum. 
These states belong to the so-called hidden spectrum [l[ . 
To find them, it is convenient to look for solutions of the 
eigenvalue equation ([3]) of the form 
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corresponding, respectively, to the eigenvalues of the 
lower hidden spectrum cr_ = — 2t coshry and to the eigen- 
values of the upper hidden spectrum er + = 2t cosh r\. In 
Eqs. A and B are normalization constants and 

X] = l/£ is a parameter accounting for the localization of 
the states. By rewriting Eq. (JTUJ) for h = irj, one has to 
solve: 

(p — 2) cosh(7y) — pcoth [i](L + 1)] sinh (77) + 2 cos(fco) = 0, 

(12) 

together with the condition that J2 X \ipo(x, y)\ 2 = 1. 
In the thermodynamic limit, L — * 00, Eq. (|12[) becomes: 
(p — 2) cosh(ry)— psinh(?7) + 2cos(fco) = 0, which is solved 
by 77 = log [p — 1 + 2 cos fco(p — 1 + cos 2 fco) + 2 cos 2 fco] , 
yielding, for the lower hidden spectrum: cr_ = 

-t (p - 2) cos(fco) + p^Jp - 1 + cos 2 (fc ) /(p-1). The 

lowest energy level is obtained for cos(fco) = 1; one has 
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Solving the eigenvalue equation ((3J for E = Eq, one ob- 
tains the wavefunction of the localized ground-state: 
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£ = l/log(- v /p+ 1) provides an estimate of the ground- 
state localization. The ground state is localized along the 
backbone and it decreases exponentially along the arms. 
When p = 2, one immediately recovers all the known 
results of the comb graph [J, y]. For each value of fc 
(cos(fco) > 0) there is a solution of (TT2]) with a differ- 
ent energy in the interval [Eq, —2t\. In a finite star-comb 
with N s (pL + 1) sites there are N s solution of this type 
and for N s — > 00 these solutions fill densely the inter- 
val [Eo,—2t]. Analogously, when cos(fco) < there are 
N s solution of (fl2|) corresponding to the spectral region 
at high energy (E = [2t, \E \]), where the other hidden 
states appear. In Fig. we plot the ground-state wave- 
function for p = 10 for a single star (dashed line) and for 
a star-comb (solid line) as a function of the distance from 
the backbone. Figure^ evidences that adding stars en- 
hances the localization of the wavefunction around the 
center of each star. 

The thermodynamic properties of non-interacting 
bosons hopping on a star-comb graph evidence also in 
this case a topology induced spatial BEC even if d < 2 
U • To elucidate this phenomenon, it is most conve- 
nient to introduce the macrocanonical ensemble to de- 
termine the fugacity z = e^ M-£;o ' as a function of the 
temperature [l[: the equation determining z is given by 
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In Eq. (fT5"]) , d(E) is the degeneracy of each single-particle 
eigenstate of the Hamiltonian (fl]) and (3 — 1/fegT. The 
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FIG. 2: (a) The normalized single-particle ground-state wave- 
function for bosons hopping on a star and star-comb graphs 
as a function of the distance y from the center. The number 
of arms is p = 10. (b) Critical temperature T c (in units of 
ks/Ej) as a function of the number of arms p. The solid lines 
are for the star-comb graph, whereas the dashed lines are for 
the single star. 



sum in Eq. (1 1 5[) is over the entire spectrum a. For free 
bosons hopping on a star-comb graph, one has 
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where N Eo (L,T), N a _ (L, T) and N a+ (L,T) denote, re- 
spectively, the number of particles at a certain temper- 
ature T in the ground-state and in the two regions <r_ 
and cr+ of the hidden spectrum. p(E), with E e uq, 
is the energy density of states of the linear chain, i.e., 
p = l/(7T-\/4i 2 — E 2 ), with (To being the region of the 
spectrum corresponding to delocalized states. The last 
term of the right-hand side of Eq. (fT6|) is the number of 
bosons in the delocalized (chain- like) states. The pres- 
ence of the hidden spectrum changes the behavior of the 
integral evaluated in the interval {—2t,2t}, since it re- 
duces it to the one describing non-interacting bosons on 
a linear chain with an impurity in one of the sites. As a 
result, letting z — > 1, the integral converges even at finite 
temperatures making possible a spatial BEC in d < 2. 

If one defines T c as the critical temperature at which 
BEC occurs, then for any T < T c the ground-state is 
macroscopically filled. Since, at the critical temperature, 
Ne (T c ) = N a± (T c ) = 0, the equation allowing to deter- 
mine T c as a function of the parameters / and t reads: 
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Equation (|17[) can be solved numerically for any value 
of /. When / ^> 1, one may expand the exponential in 
Eq. (fT7|) to the first order in the inverse of the critical 
temperature T c . Upon inserting the value of the ground- 
state energy (IT31 in Eq. (fTT)) . one easily finds that the 



critical temperature T c is given by 

tf Vp( 2 + Vp) 



(18) 



Equation (|18|) has been checked numerically and it is in 
excellent agreement with the numerical solution of Eq. 
(fT7| for / 3> 1, the error being of order 1/f. If one 
compares the critical temperature of the single star and 
of the star-comb graphs, one immediately sees that, once 
the number of arms is fixed, T c is always enhanced in the 
latter realization. This is shown in Fig. [2]d. 

One may now use Eq. (|18[) to determine the con- 
densate fraction as a function of the scaled temperature 
T/T c . In the thermodynamic limit, the number of parti- 
cles in the delocalized states is given by 
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In the last equation the exponential has been expanded 
to the first order in (3: this approximation holds for / 3> 
1 and it is in very good agreement with the numerical 
evaluation of the integral in Eq. (fH?|) also in a large 
neighborhood below T c . From Eqs. (HHJ) and (fT§|) one 
gets the number of particles in the localized states No = 
N Eo + N a _: the fraction of condensate, for T < T c , is 
then given by 
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For / ranging from 10 3 to 10 9 , the results provided by Eq. 
(|20f differ from those obtained by the numerical evalu- 
ation of Nq from Eq. p^|) by less than 1%. Equation 
([2TJ)l clearly shows that the condensate has dimension 1 
just as cigar-shaped one-dimensional atomic Bose con- 
densates @- 

Due to the topology induced spatial condensation on 
the backbone, one should expect an inhomogeneous dis- 
tribution of the bosons along the arms of the network. 
The average number of bosons Nb at a site (x, y, z) de- 
pends - due to the symmetry of the graph - only on the 
distance y from the backbone. It is not difficult to show 
(see also Ref. Q) that, away from the backbone (y 1) 
and once the filling is fixed, Nb depends only on the 
scaled temperature T/T c and it is given by 



N B {y;T/T c ) _ T 
f ~ T c 



(21) 



Topology induced spatial BEC in a system of non- 
interacting bosons hopping on a star-comb graph predicts 
then a rather sharp decrease of the number of bosons at 
sites located away from the backbone. The linear de- 
pendence is consistent with the observation that, in this 
system, the condensate has dimension 1. 
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In this paper we have investigated in detail another sit- 
uation where the role played by the network's topology is 
crucial in determining the thermodynamic properties of 
the system. We focused on the properties of free bosons 



hopping on a star-comb network, finding an enhancement 
- with respect to other inhomogeneous structures previ- 
ously investigated [J, H, 0, 0] - of the critical temperature 
at which the particles condense. 
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